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ABSTRACT 

We calculate simulated images of disks perturbed by embedded small planets. 
These 10—50 M bodies represent the growing cores of giant planets. We examine 
scattered light and thermal emission from these disks over a range of wavelengths, 
taking into account the wavelength-dependent opacity of dust in the disk. We 
also examine the effect of inclination on the observed perturbations. We find 
that the perturbations are best observed in the visible to mid-infrared. Scattered 
light images reflect shadows produced at the surface of perturbed disks, while 
the infrared images follow thermal emission from the surface of the disk, show- 
ing cooled/heated material in the shadowed/brightened regions. At still longer 
wavelengths in the sub-millimeter, the perturbation fades as the disk becomes 
optically thin and surface features become overwhelmed by emission closer to- 
ward the midplane of the disk. With the construction of telescopes such as TMT, 
GMT and ALMA due in the next decade, there is a real possibility of observ- 
ing planets forming in disks in the optical and sub-millimeter. However, having 
the angular resolution to observe the features in the mid-infrared will remain a 
challenge. 

Subject headings: planetary systems: formation — planetary systems: proto- 
planetary disks — radiative transfer 



1. Introduction 



Giant planet formation remains an unsolved problem. The two leading theories are disk 
instability and core accretion. In disk instability, the protoplanetary disk gravitationally 
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fragments into Jupiter-mass clumps which condense into planets. In core accretion, smaller 
bodies agglomerate into larger and larger bodies until they are massive enough to accrete a 
gaseous envelope. There are theoretical arguments for and against both scenarios, but the 
aim here is to identify signatures of planet formation in order to provide direct observational 
evidence for one scenario or the other. Simulated scatte red light signatures of plane t forma- 
tion by disk instability have been discussed elsewhere (IJang-Condell fc Bossl 120071 ). so the 
focus of this paper is on signatures of core accretion. 

Giant planets forming by core acc retion need to have co res of 10-20 M e to be massive 
enough to accrete a gaseous envelope (IHubickyj et al.l 120051 ). For reference, Jupiter is ap- 
proximately 300 M e . This predicts that sizeable planet embryos form before circumstellar 
gas disks dissipate. These disks are typically modeled as passive accretion disks, meaning 
that the dominant heating source is stellar irradiaion, although there still is sufficient viscous 
dissipation for the disk to continue accretion at a slow rate onto the star (< lO _7 M yr _1 ) 
and for a modest amount of heating at the disk midplane. These disks are optically thick 
and gas-dominated. Their temperature structure is strongly dependent on heating from 
stellar irradiation, particularly on the angle o f incidence of the star l ight on the disk surface 
dChiang fc Goldreichl fl997l : ICalvet et al.lll99ll : b'Alessio et~aDll998l . Il999l . boOlh . 



We examine planets forming at distances of 1—8 AU from the stars. Interior to 1 AU, the 
innermost regions of protoplanetary disks are expected to become opticall y thin because of 
the sublimation of dust grains, creating a hot inn er rim at ^ 0.1 — 0.5 AU (IDullemond et al. 



200 ll ; llsella fc Nattal 120051 : iD'Alessio et al.l 120061). This is supported by observ ations of T 
Tauri stars in the near- to mid-IR (jEisner et al.ll2009l . l2005l : lMuzerolle et al.ll2003l ). Although 
larger radii for inner holes in protoplanetary disks are observed, this is typically for higher 
mass Herbig Ae/Be stars, which are both hotter and more luminous than T Tauri stars, 
resulting in the evaporation of dust grains at larger distances from the star. The typical 
inner disk radius for T Tauri stars is below 0.5 AU, which is the minimum disk distance 
modeled in this study. For this reason, the emission from the inner disk rim can be treated 
as a separate problem from the study presented in this paper. 

Other work on simulated images of growing planets focuses on Jovian-mass planets 
rapidly accreting material from their surroundings disks. These planets are large enough to 
create large-scale structures in disks such as gaps or spiral densit y waves, and the planets 
them selves are observable because they are heated by accretion (e.g. IWolJl2008l : lKlahr fc Kiev 
20061 ). Planet formation has already completed and it is difficult to determine the mecha- 
nism by which it occured. By contrast, the study presented in this paper focuses on planets 
which are too small to open gaps or accrete significant amounts of gas (10 — 50 M e ). 



The planet-disk models adopted here are presented in Uang-Condelll (120081 ) and sum- 
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marized in this section. The planet is predicted to gravitationally compress the disk in the 
vertical direction, creating a shadow paired with a bright spot, leading to temperature vari- 
ations. We calculate both scattered light and thermal emission from the perturbed disk for 
planets of mass 10 — 50 M e at distances of 1 — 8 AU from their host stars. 

In §21 we describe the procedure we use for modeling radiative transfer in the disk. In 
§3] we show the results of the modeling for the full suite of planet-disk models over a range 
of wavelengths. In §0] we discuss the potential observability of the planet perturbations we 
have modeled. In §5] we present our conclusions. 

2. Model Description 
2.1. Disk Model: Thermal Structure 



The mo del for therma l pertu rbations in the vicinity of an embedded protoplanet are 
presented in Ijang-Condelll (120081 ). The stellar parameters are mass M* = 1M , radius 
R* = 2.6 i? , and effective temperatur e = 4280 K. The radius and effective temperature 
correspond to a stellar age of 1 Myr (jSiess et al.ll2000h . We assume a viscosity parameter 
of a = 0.01 and accretion rate M = lO _8 M yr _1 , which are typical of T Tauri stars 



flHartmann et al.lll998l : iD'Alessio et al.lll998h . 



The model for the disk is a based on an accretion disk model with viscosity parametrized 
by a dimensionless constant a, so that the viscosity is v — ac s H where c s = y/ kT / fh is 
the the rmal sound speed at the midplane of the disk and H = c s /Q is the pressure scale 
height ( IShakura fc Sunyaevl Il973l ). Here, k is the Boltzmann constant, T is the midplane 
temperature, fh is the mean molecular weight of molecular hydrogen, and Q is the Keplerian 
angular velocity. For the low accretion rate adopted and the distances in the disk in which 
we are interested, stellar irradiation heating dominates, though we include viscous heating 
for completeness. 

The amount of heating from stellar irradiation at the surface depends on the angle of 
incidence of the stellar rays. Based on an analytic solution to the 1-D plane-parallel case, we 
calculate surface heating semi-analytically in a 1+2D manner, where we numerically integrate 
heating contributions over the shape of the disk surface. We include the differential rotation 
of the disk by calculating the rate of radiative cooling and heating as disk material enters or 
exits the shadowed or brightened regions. 



We assume vertical hydrostatic equilibrium in the disk, including contributions to the 
vertical gravity from both the star and the planet. We treat the planet embedded in the disk 
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as a point mass, whose gravitational potential compresses the disk in the vertical direction, 
creating a dimple at the surface. The dimple at the surface results in a shadow at the planet's 
position relative to the unperturbed disk, while the outer edge is brightened. 

The radiative transfer calculations rely on a separation of the radiation into two regimes: 
the short wavelength spectrum of the stellar emission, and the long wavelength spect r um o f 



the thermal emission of the disk material. We use the opacities from iD'Alessio et al.l (120011 ) 
using a dust model with parameters a max = 1 mm, T = 300 K, and p = 3.5, assuming that 
the dust opacities are constant throughout the disk. The values for the opacities (in cm 2 g -1 ) 
are as follows: the Rosseland mean opacity is xr — 1-91, the Planck mean opacity integrated 
over the disk spectrum (300 K) is ftp = 0.992, and the Planck mean opacities integrated 
over the stellar spectrum (4000 K) are n* P = 1.31 for absorption alone and Xp — 5.86 for 
absorption plus scattering. The absorption fraction is then = ftp/xp, while the scattered 
fraction is a = 1 — a^. The Rosseland mean opacity is used to calculate the photosphere 
of the disk, and Xp 18 use d to calculate the surface of the disk. We assume that the dust is 
well-mixed with the gas and constant throughout the disk. 

The initial conditions vary slightly in the overlap regions between simulations as the 
size and position of the simulation box vary. The sizes of the boxes are scaled to the Hill 
radius, 

f Tfi \ 

rHm = UtJ a - (1) 

Thus, the boundary conditions depend on the mass and the distance of the planet. In order 
to make a fair comparison between the various planet models, we assemble a composite disk 
from the initial conditions of the largest boxes at each distance from the star. We interpolate 
the densities and temperatures smoothly across the overlap regions. The magnitude of the 
variations in temperature and density should roughly scale, so we normalize the temperatures 
and densities of each planet-disk model to this composite set of initial conditions. 



2.2. Opacities 



In order to calculate simulated images of the disk, we need a good model for the opacities. 
Although mean opacities were good enough for calculating the thermal structure of the 
disk, to calculate images of the disk, we need the frequency-dependent absorption and total 
extinction coefficients for the dust model. 

The opacities were calculated using a Mie scatte ring code develope d bv lPollack & Cuzzi 



(Il980l ). The composition of the dust is that used in lPollack et al.l (11994 ) . consisting of water, 
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troilite, astronomical silicates, and organics. We adopt a size distribution for the dust of 
of n(a) oc a -3 5 where a is the radius of the grain, with maximum and minimum radii of 1 
mm and 0.005 microns, respectively. The large maximum grain size represents coagulation 
of grains in protoplanetary disks. 

In FigurefU we show the absorption (dotted line) and total extinction (solid line) coeffi- 
cients in units of cm 2 g _1 versus wavelength for our adopted dust model. The total extinction 
includes both scattering and absorption. If k v and \v are the frequency-dependent absorp- 
tion and extinction coefficients, respectively, then the albedo is ratio of scattering to total 
extinction, so uj v = 1 — K v jxv Since we assume that the dust is well-mixed, the dust density 
simply scales with the gas density. 



2.3. Scattered Light Imaging 

We calculate scattered light from the disk assuming that it is face on and that stel- 
lar irradiation is primarily scattered from part of the disk that is optically thin to stellar 
irradiation. The optical depth to stellar light at a frequency v is 



fxupdl (2) 
Jt 



where \v 18 the extinction coefficient, p is the density of the disk, and t is the path between 
the star and the point in the disk. From this point on, the dependence on frequency, z/, will 
be implicit and we will omit the subscripts. We define the surface of the disk, z s {r) to be 
where t* = 2/3. The cosine of the angle of incidence at radius r is 

a(r) = & - Zs,r) + 4it> * ( 3 ) 

[{z f s f + 1] 1/2 (1 + z 2 Jr 2 ) 1/2 3tt (r 2 + z 2 ) 1/2 

The first term is a purely geometric term, and the second term accounts for the finite size 
of the star, in constrast to illuination by a point source. The slope of the surface is defined 
as z' s = dz s /dr. We use a log-log interpolation to calculate both z s and z' s as a function of r. 

Assuming isotropic scattering, the brightness of dis k as a function o f frequ ency, I*, can 



be calculated by integrating the following equation from D'Alessio et al. ( 19991 ): 

dP ( X - k)RI 



dZ 4/ 2 



-B(T*) exp[-r* - r ohs (Z)} (4) 



where I is the distance to the star, B(T*) is the Planck function at the effective temperature 
of the star, and r D bs is the optical depth to the observer, r b s = J obs Xupdl, integrated along 
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the line-of-sight from the observer to the disk. For a face-on disk, this is integrated along 
the vertical axis, so Z = z. 

At low optical depths, p is small, so contributions to I s are small. At high optical 
depths, the exponential term goes to 0, so those contributions are also small. Only the 
regions close to the surface contribute significantly to the scattered light intensity, so we 
assume I w \Jr 2 + z 2 s . We assume that the disk is locally plane parallel at the surface, so 
that if % is the angle between the line of sight to the observer and normal to the surface, 
t* = coszT bs(^)//i. Note that % is not the inclination of the disk, but rather the angle with 
respect to the disk surface, which is not generally parallel to the ecliptic. Then Eq. 
becomes 

^Tz =d i = ipf^) B(r * )exp[ - (1+cosWr] - (5) 

Integrating from r = — > oc, 



rscatt _ upR*B(T*) 



4(r 2 + zl)(p + cosi) 



(6) 



2.4. Thermal Emission 

At longer wavelengths, thermal emission from th e disk material become important. The 
thermal emission from a disk is calculated following iD'Alessio et al. by integrating 

the equation 

dP 

— = KpB(T d ) exp(-r obs ) (7) 

along the line of sight, where is the local disk temperature and B(Td) is the thermal 
emission at that temperature. The thermal emission and scattered light are summed to give 
the total disk surface brightness: 

/ disk = /S + / t > ^ 



3. Results 

In Figures l2lfTT)l we show the results of simulating images planets in disks from 1 to 300 
microns. Face-on disks (0° inclination) with planets are displayed in Figure [3 Each disk 
image shows a different wavelength, with 10, 20 and 50 M planets embedded at 1, 2, 4 and 
8 AU. The 10 M planets are at the 9 o'clock position and the 20 M planets are are the 12 
o'clock position. The 50 M e planets are in the lower right quadrant of the disk. They are 
offset from each other because their simulation boxes overlap slightly in the radial direction. 
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In subsequent figures, we allow the boxes to overlap because the perturbation inward of the 
planet is minimal. 

The shadows are roughly centered at the planet's position in radius and azimuth, and 
brightened spots appear on the far side of the dimple. We define the contrast to be c = 
f/fo — 1 where / is the radially-dependent unperturbed disk brightness and / is the disk 
brightness in the presence of the planet. The depths of the shadows /bright spots regions in 
the face-on disk can then be characterized by the minimum/maximum value of the contrast. 
The sizes of the shadows/ bright spots are characterized by the area enclosed by a contour 
traced at half the minimum/maximum constrast. These values are tabulated in columns 5, 
6, 8 & 9 of Tabled] by planet mass (column 1), orbital distance (column 2), and wavelength 
of observation (column 4). The "spot radius" is the equivalent radius of the half-min/max 
areas A, 



In Figure [3l we plot the maximum/minimum values of the contrast versus wavelength, 
with line type and color indicating orbital distance and symbols indicating planet mass. The 
contrasts generally increase with mass, while the dependence on distance in more compli- 
cated. The wavelength dependence will be discussed in further detail below. 

Figure [H shows the dependence of r spot /rHiii with wavelength, with symbol shape and 
color indicating mass, and horizontal offset indicating distance. The clustering of the points 
indicates that r spot scales roughly with rnm, with brightened spots being generally larger 
than shadows. In other words, the size of the shadow scales roughly with the distance from 
the star and the cube root of the planet mass. 

In Figures El El IHl l9landfT0lwe show the effect of inclincation on the appearance of the 
perturbed disk at wavelengths of 1, 3, 10, 30, 100, and 300 microns, respectively. Each figure 
shows, from left to right, an angle of inclination of 30°, 45° and 60°. From top to bottom, 
the planet mass is 10, 20, and 50 M®. The inclination axis is always pointing upwards, so 
that the upper part of the disk is further from the observer than the bottom part. In each 
disk model, planets are inserted at 1, 2, 4, and 8 AU, along each minor and major axis to 
demonstrate how the phase of the planet affects the appearance of the perturbation. 




(9) 



3.1. Scattered light 



In the optical to near-IR (1 and 3 fxm, upper plots in Figure ED, the image of the 
disk is predominantly from scattered light. These images trace the contours of the surface 
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of the disk. The shadowing and illumination are clearly visible. For reference, the stellar 
brightness is 1.54(rf/100 pc) -2 Jy at 1 micron, and 0.77(d/100 pc) -2 Jy at 3 microns, where d 
is the distance to the star. At these wavelengths, the brightness of the star overwhelms any 
emission from the disk, so very good starlight suppression, i.e. coronography, is necessary to 
resolve any features in the disk. Coronography is also a good way to suppress emission from 
hot inner walls of the disk, since this may be a significant source of near-IR flux in the inner 
0.5 AU of the disk. 

The shadow contrast is generally deeper than the bright contrast because it is easier 
to create a depression than it is to puff up disk material due to heating, espec i ally i n the 



surface layers of the disk where the density is very low. As seen in lJang-Condelll ( 120081 ). the 
heating from illumination is greater interior to 8 AU, because the densities are higher in the 
inner regions. The exception to this is for the largest planet at 2-4 AU, where the bright 
contrast rises above 1, whereas the shadow contrast has a floor at —1. On the other hand, 
the area of the brightened regions are generally larger than the shadowed region, except for 
the smallest planet at 8 AU, where the amount of brightening is very small. 

The effect of inclination on planet dimples can be seen in Figures [5] and [61 One general 
effect is that as the angle of inclination increases, the brightness of the disk also increases. 
Examining equation ([6]), we see that as z, the angle of observation with respect to the surface 
normal, increases, cosi decreases, increasing the overall value of I s . The interpretation for 
this is that the observer is viewing more scatterers along the line of sight through the disk, 
thus the disk appears brighter. Since the disk images are tipped up so that the bottom half 
is closer to the observer, the value of i is greater, so the bottom halves appear brighter. The 
planet perturbation on top of this general effect causes the rebrightened region of the disk 
to appear brighter still when it is on the near side. 



3.2. Mid-IR: thermal emission from the surface 

The mid-IR (10 and 30 /xm, center plots in Figure [2]) probes thermal emission from 
the surface layers of the disk. At these wavelengths, the observed features are caused by 
temperature perturbations resulting from cooling in the shadowed regions and heating in 
the illuminated regions. The size and contrast of the perturbations show the same variation 
with planet size and distance as seen in the scattered light images. Again, the shadows 
appear deeper than the brightened spots, except when the bright contrast > 1.8 because of 
the natural floor in the shadow contrast at —1. The stellar brightness is 0.107(rf/100 pc) -2 
Jy at 10 microns, and 1.34 x 10 _2 (d/100 pc) -2 Jy at 30 microns. 
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The greatest amount of contrast occurs at 10 /im. The 10 fim image also shows the 
most radial contrast in the disk overall because the outer regions of the disk are too cold to 
emit effectively at this wavelength: 10 /im is shortward of the exponential cutoff in thermal 
emission in the outer disk. The steep background radial gradient of the unperturbed disk 
itself will make it difficult to isolate perturbations from embedded planets. 

In the inclined disks (Figures [7] and [8]) there appears to be a brightening associated 
with higher inclinations, similar to the scattered light images, but of lesser magnitude. This 
brightening is caused by an optical depth effect, because as more heated or cooled material 
aligns along the line of sight, the stronger or weaker the thermal emission becomes. 



3.3. Far-IR to radio: thermal emission from the midplane 

At sub-millimeter wavelengths, the stellar brightness is negligible compared to the ther- 
mal emission from the disk. In this regime the perturbations induced by the embedded 
planets are hardly visible (0.1 mm and 0.3 mm, lower plots in Figure [2]). The contrast is 
modest compared to shorter wavelengths. The scatter in r spot at 300 /im in Figure H] is largely 
because the contrast is so modest that the half-min/max contours are hard to constrain. 

The low contrast arises because the thermal structure of the disk is affected primarily 
in the upper layers of the disk. As the disk becomes optically thick to its own thermal 
radiation, the effects of small perturbations at the surface become washed out. The contrast 
between the inner and outer disk also becomes less pronounced because these wavelengths 
are in the Rayleigh- Jeans regime of the blackbody spectra. 

However, other effects not modeled in this study may very well create observable signa- 
tures. In particular, dynamics within the Hill sphere is not properly included. Circumplane- 
tary disk formation and accretion onto the forming planet will create a density enhancement 
and heating at the midplane, which will increase the thermal emission in this region, particu- 
larly at wavelengths where the disk becomes optically thin. Better models need to be carried 
out to properly model planet signatures in the radio, particularly to determine whether small 
planets create shadows versus bright enings. 



4. Discussion: Observability 



In order to observe the phenomena modeled in this paper, we need to look to nearby 
star-forming regions to find sufficiently young pr otoplanetary disks. T he Taurus- Auriga star- 
forming is one of the nearest, at about 140 pc (IKenyon et al.lll994l ). At this distance, a 1 
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AU perturbation subtends 0.007". 

The best wavelengths for observing planet shadows as modeled in this paper are in the 
visible to mid-infrared. Scattered light images are sensitive to small perturbations in density 
at the surface of the disk. However, this is both an advantage and disadvantage. Small 
variations such as those caused by small planets can cause large shadowing effects. On the 
other hand, density variations might be caused by other processes, such as outflows, disk 
turbulence, and dust dumpiness, just to name a few. Differentiating planet-induced shadows 
from those caused by other processes may prove to be difficult. 

At 1 /xm, a diffraction limit of X/D = 7 mas requires a telescope diameter of 29 me- 
ters. This is larger than any current ground-based telescope or planned space telescope. 
However, with the development of large segmented-mirror telescopes such as the Giant Mag- 
ellan Telescope (GMT)0 and the Thirty Meter Telescope (tmt||, these surface shadows 
may be resolvable in the near future. However, scattered light observations of the planet 
perturbations requires very good suppression of the starlight at small inner working angles. 

For mas resolution at longer wavelengths, single-dish telescopes are increasingly infea- 
sible, but interferometry is a promising alternative. The Center of High Angular Resolution 
Astronomy (CHARA) Arrajl^ at Mt. Wilson Observatory has baselines of up to 330 meters 
and operates in the optical and 2.0 — 2.5 jim range, allowing for sub-mas resolution. The Very 
Large Telescope Interferometer (VLTljfl can achieve up to 2 mas resolution in J, H and K 
bands (1 — 2.4 fim) with the AMBER instrument. VLTI-MIDI is a two-way beam combiner 
at N-band (8 — 13 /im), for an angular resolution of about 10 mas. A next-generation mid- 
IR instrument for VLTI, MATISSE0, is also under development for L, N, M, and Q bands 
(3 — 25 fim). The challenge for these instruments is covering sufficient u-v space to build up 
an image, achieving sufficient dynamic resolution so that the star does not overwhelm the 
disk emission, and having enough sensitivity to detect the planet shadows. 

Mid-infrared wavelengths, which probe temperature perturbations caused by shadowing 
and illumination at the surface, are better for identifying planet perturbations. This is 
because shadows and brightenings caused by transitory phenomena will not persist as long 
those caused by an embedded planet, and will not result in strong cooling and heating in 



1 http://www.gmto.org/ 
2 http:/ /www. tmt.org/ 

3 http:/ /www. chara.gsu.edu/CHARA/index.php 

4 http : / / www.eso .org / sci / facilities / par anal / telescopes / vlti / index, html 
5 http:/ /www.oca.eu/matisse/ 
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the surface layers. Although the shadow contrast is good at 10 /im, the dynamic range 
between the inner and outer disk is quite high, over 4 orders of magnitude from 0.5 to 10 
AU. Dynamic range is less of an issue at 30 fim. The 30 /xm images shown in Figures [2] and 
[8] are very similar to the 1 /xm images in Figures [2] and [5j these two different wavelengths 
trace the planet perturbations equally well. 

A resolving power of 7 mas at 30 /xm requires firstly getting above the Earth's atmo- 
sphere, then having a baseline of 870 m: this would almost certainly require a formation- 
flying space mission. However, there are no near-term prospects for such a mission. 

In the radio, the Atacama Large Millimeter/submillimeter Array (ALMAfl will have 
mas resolution when it reaches its full operating capacity in 2012. While the resolution 
would be good enough to resolve a planet shadow, the perturbations have very low contrast 
at 0.3 mm, the shortest wavelength that ALMA will operate at. The problem here is one of 
sensitivity rather than angular resolution. 

This study was restricted to planets forming interior to 10 AU, based on scenarios of 
giant planet formation within our own Solar System. However, discoveries of planets around 
HR 8 799 and Fomalhaut i ndicate that exopl anets can form tens to hundreds of AU from their 



stars (IMarois et al.ll2008l : iKalas et al.ll2008l ). The size of planet shadows in protoplanetary 
disks scale roughly with the planet-star distance. So there is hope that planets forming at 
these far distances from their stars might be observable in the near future, perhaps even 
with current instrumentation. 



5. Conclusions 

Observable perturbations from planets roughly scale with planet mass and orbital dis- 
tance. The features are most easily observed in visible to mid-infrared wavelengths. The 
visible regime probes scattered light off the surface of the disk, revealing shadows caused by 
the gravitational potential of the planet. Mid-IR wavelengths reveal cooling/heating from 
shadowed/brightened regions at the surface. Confirmation of shadows and brightenings seen 
in scattered light images of disks should therefore be carried out at near- to mid-IR wave- 
lengths. In longer wavelengths toward the sub-millimeter, the contrast of the perturbations 
becomes smaller as the disk becomes optically thin. 

The validity of the simulations presented in this study is limited by the fact that hy- 
drodynamics are not included. This means that large-scale effects such as gap-clearing are 



6 http: / / almaobservatory.org 
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not well-modeled. If an annular gap does form, this density perturbation will be subjected 
to shadowing and illumination in the same way local dimples are, so radiative feedback on 
gaps should enhance these structures as well. Thus, the local planet shadows modeled in this 
paper represent lower limits on perturbations to disks that might be caused by embedded 
plaent cores. 

Another missing physical process is accretion onto the forming planet core. This process 
is difficult to capture even in high resolution three-dimensional hydrodynamic simulations 
because of the extremely high dynamic range necessary to resolve the final capture of gas into 
the atmosphere of the growing planet. Heating from accretio n onto the planet could alter 



the results of this study. Simulations by lHubickyj et al.l (120051 ) indicate that growing planets 



accrete slowly over a long period of time, but once it crosses a threshold mass of roughly 
10 — 20 M®, it accretes gas rapidly over a short time period. Thus, accretional heating is 
likely to be only transitionally important. On the other hand, planets might grow so rapidly 
that 50 Earth mass planets might never be seen, being only a transitional stage on the way 
to becoming a Jovian- type planet. 

The author thanks John Debes for helpful discussions and advice in developing the dust 
opacity model used in this paper; Lee Mundy and an anonymous referee for helpful comments; 
and Matthew Condell for assistance with software development. This work was performed 
under contract with the Jet Propulsion Laboratory (JPL) funded by NASA through the 
Michelson Fellowship Program. JPL is managed for NASA by the California Institute of 
Technology. 
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Fig. 1. — Wavelength-dependent opacities for the dust model adopted for calculating sim- 
ulated images. Dotted line: absorption only. Solid line: the total extinction (absorp- 
tion+scattering) . 
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Fig. 2. — Planets of varying masses embedded in a face-on disk at 1, 2, 4 and 8 AU seen at 
varying wavelengths. From left to right, top to bottom: 1, 3, 10, 30, 100, and 300 microns. 
Planets at the 9 o'clock position are 10 M e , planets at the 12 o'clock position are 20 M e , 
and planets at the 3-5 o'clock positions are 50 M e . 
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A(microns) 

Fig. 3. — Minimum/Maximum contrast versus wavelength. Orbital distances of 1, 2, 4, and 8 
AU are represented by black dotted, green short- dashed, blue long-dashed, and solid magenta 
lines, respectively. Planet masses are indicated by crosses/squares/stars for 10/20/50 M®. 
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Fig. 4. — Equivalent radii r spot of shadowed/brightened regions, scaled to the Hill radius rnm? 
versus wavelength. The points have been horizontally shifted according to to distinguish 
between orbital distances, but the wavlengths correspond to 1, 3, 10, 30, 100, and 300 
fim exactly. Black triangles/green squares /magenta diamonds indicate planet masses of 
10/20/50 M e . Filled symbols indicate shadows and open symbols indicate bright spots. 



- 18 - 




AU AU AU 



Fig. 5. — Simulated images of inclined disks perturbed by embedded planets at 1 fim. The 
stellar brightness at 1 /im is 1.54 Jy at a distance of 100 pc. Each image shows how a planet 
of a given mass will perturb the appearance of the disk at 1, 2, 4, and 8 AU, and at phases 
0, 7r/2, 7r, and 37T.2. From left to right, the disks are inclined at 30, 45, and 60°. The planet 
masses from top to bottom are 10, 20 and 50 M@. 




Fig. 6. — Simulated images of inclined disks perturbed by embedded planets at 3 fim. The 
stellar brightness at 3 fim is 0.77 Jy at a distance of 100 pc. See Figure [5] for description. 
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Fig. 7. — Simulated images of inclined disks perturbed by embedded planets at 10 fim. The 
stellar brightness at 10 fim is 0.107 Jy at a distance of 100 pc. See Figure [5] for description. 




Fig. 8. — Simulated images of inclined disks perturbed by embedded planets at 30 fim. See 
Figure [5] for description. 




Fig. 9. — Simulated images of inclined disks perturbed by embedded planets at 0.1 mm. See 
Figure [5] for description. 




Fig. 10. — Simulated images of inclined disks perturbed by embedded planets at 0.3 mm. 
See Figure [5] for description. 



-24- 



Table 1. Contrast and size of shadowed and brightened regions. 
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Table 1 — Continued 
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